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Abstract 

A general method for the derivation of asymptotic nonlinear shallow water 
and deep water models is presented. Starting from a general dimensionless 
version of the water-wave equations, we reduce the problem to a system of two 
equations on the surface elevation and the velocity potential at the free surface. 
These equations involve a Dirichlet-Neumann operator and we show that all 
the asymptotic models can be recovered by a simple asymptotic expansion of 
this operator, in function of the shallowness parameter (shallow water limit) 
or the steepness parameter (deep water limit). Based on this method, a 
new two-dimensional fully dispersive model for small wave steepness is also 
derived, which extends to uneven bottom the approach developed by Matsuno 
|27| and Choi [28]. This model is still valid in shallow water but with less 
precision than what can be achieved with Green-Naghdi model, when fully 
nonlinear waves are considered. The combination, or the coupling, of the new 
fully dispersive equations with the fully nonlinear shallow water Green-Naghdi 
equations represents a relevant model for describing ocean wave propagation 
from deep to shallow waters. 
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1 Introduction 



The propagation of surface waves through an incompressible homogenous inviscid 
fluid is described by the 3D Euler equations combined with nonlinear boundary 
conditions at the free-surface and at the bottom. This problem is extremely difficult 
to solve, in particular because the moving surface boundary is part of the solution. 
The complexity of this problem led physicists, oceanographers and mathematicians 
to derive simpler sets of equations in some specific physical regimes. Equations thus 
obtained may be divided into two groups, namely shallow-water models and deep- 
water models. One of the goals of this paper is to clarify the range of validity of 
these models and to show why at least two different asymptotic models are necessary 
for a correct description of ocean waves. 

In shallow water conditions, the classical approach is based on a perturbation 
method with respect to a small parameter \i = (h /L) 2 (/z the characteristic water 
depth, L the characteristic horizontal scale), in order to reduce the three-dimensional 
equation system to a two-dimensional one. This method, initially introduced by 
Boussinesq [I], allows to derive several shallow- water equations, which are named 
"Boussinesq-type equations". A large class of such equations can be expressed in 
the following 2D nondimensional form: 

r d t ( + v-(hv) = o(^) 

\ d t V + e(V ■ V)V + VC = fiD + 0{^ 2 ) K ' 

where e = a/h (a the order of free surface amplitude) is the nonlinearity parameter 
and the dimensionless flow variables are: h the water depth, £ the surface elevation 
and V the depth-averaged velocity. T> characterizes non-hydrostatic and dispersive 
effects and is a function of wave variables and their derivatives. Higher-order Boussi- 
nesq equations can be derived (e.g. Madsen et Schaffer [2]), but in this paper we 
restrict our analysis to 0(fi 2 ), which is a good approximation for most of nearshore 
wave applications (e.g. Madsen et al. [3], [1], Cienfuegos et al. [5]). 
In its classical form, Boussinesq wave theory is a one-dimensional approach based 
on the assumptions of weak dispersion, weak non-linearity and balance between dis- 
persion and non-linearity: e = O(fi) <C 1. The KdV [6] and BBM [7] equations 
can also be derived in the case of unidirectional waves; this approach is relevant to 
study fundamental wave dynamics problems, such as ID solitary wave propagation 
on flat bottom (e.g. Boussinesq [T] and Benjamin [S]). However, the classical Boussi- 
nesq assumptions may severely restrict applicability to real word wave propagation 
problems. Applications to coastal zone, have motivated theoretical developments 
(see reviews by Dingemans [9], Madsen et Schaffer [2], Kirby [10] and Barthelemy 
[TT] ) for extending the range of applicability of Boussinesq-type equations in term 
of varying bottom, dispersive and nonlinearity effects, which play an important role 
in the nearshore wave dynamics. 

The 2D Boussinesq equations (e = 0(/i) <C 1) for non flat bottom were first derived 
by Peregrine |12j : 

V = ^ V (V ■ (hd t V)) - %V 2 d t V . (2) 
2 o 

For many coastal applications the weak dispersion of these equations is a critical 
limitation. Witting [13] proposed a method based on Pade expansion of the exact 
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linear phase velocity to improve Boussinesq-type equations. From this method, 
several equations (order 0(/i 2 )) with improved dispersion characteristic have been 
derived (e.g. Madsen et al. [33], Nwogu [15], Schaffer and Madsen [16], Bona et al. 

PUDS!)- 

In 1953, a breakthrough treating nonlinearity was made by Serre (see Barthelemy 
[TT] for a review). He derived ID fully nonlinear (e = 0(1)) weakly dispersive 
equations for horizontal bottom: 

V = ^-d x {h\V xt + VV XX - (V x ) 2 )). (3) 
3h 

The same system was obtained later by Su and Gardner [19]. Seabra-Santos et al. 
[20] have provided an extension of this model to non flat bottom: 

V = -\ (d*(h 2 {\p + \Q)) + b x h{\p + Q)) , (4) 

with P = -h(V xt + VV XX - (V x ) 2 ) and Q = b x (V t + VV X + b xx V 2 ). Dingemans 
[9] expressed mistrust for the validity of these equations for uneven topographies, 
considering that the derivation required the assumption of vertical uniformity of the 
horizontal velocity. Cienfuegos et al. [21] have shown that Serre equations can be 
formally established without assumption on the velocity profile. 
Finally, Green and Naghdi [22] derived 2D fully nonlinear weakly dispersive equa- 
tions for uneven bottom which represents a two-dimensional extension of Serre equa- 
tions. Except for being formulated in terms of the velocity vector at an arbitrary 
z level, the equations of Wei et al. [23] are basically equivalent to the 2D Serre or 
Green-Naghdi equations; this is also the case of the equations derived in [23] through 
Hamilton's principle. 

The range of validity of all the models introduced above may vary as far as the 
nonlinearity parameter e is concerned, but they all require that the shallowness 
parameter fi is small. In deep water (/i <^ 1) the Boussinesq wave theory fails but 
it is yet possible to derive asymptotic expansions from the water-wave equations 
under the condition that the steepness e^JJi = a/A is small. Such an approach was 
first made in ID and flat bottoms by Matsuno [25] and extended to uneven bottoms 
[2"E] . 2D weakly transverse waves |27j and higher order expansions [28J. Since one 
always has e < 1, the small steepness assumption eJJi <C 1 is also satisfied in the 
shallow water regime /i < 1 discussed above and this is the reason why it is often 
claimed that the models derived in [251 [26| [271 [28] are valid in the whole range 
fi E (0, oo). However, as we show here, their precision is then far below the one of 
the Green-Naghdi equations. 

In this paper we propose a systematic derivation of all the models evoked in this in- 
troduction (and of some new ones); to this end, we use the global method introduced 
in the recent mathematical work |29j. Starting from a general nondimensionalized 
version of the water-wave equations (which takes into accounts the different nondi- 
mensionalizations used in deep and shallow water), following [30, El], we reduce 
the problem to a system of two equations on the surface elevation and the velocity 
potential at the free surface. These equations involve a Dirichlet-Neumann operator 
and we show that all the asymptotic models can be recovered by a simple asymp- 
totic expansion of this operator. In the deep water setting (/i ^ 1) this expansion 
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is made with respect to e (which is small, as a consequence of the small steepness 
assumption e-JJl <C 1 and the deep water assumption /i 1); in the shallow-water 
setting, we rather choose to expand the Dirichlet-Neumann operator in terms of [i 
since e is not always small (e.g., e ~ 1 for Green- Naghdi). 

As said above, all the asymptotic models are obtained in a systematic way with the 
present approach. In addition, we derive a new fully dispersive model for 3D waves 
over uneven bottom. Our method brings clarifications about the validity domain of 
the main asymptotic wave models. In particular, 

• We can easily show that the shallow water limit of the Dirichlet-Neumann 
expansion corresponding to the deep-water models has a very poor precision 
in the fully nonlinear case e ~ 1; this proves that the deep-water models such 
as those of [251 12SI 123 [28] are not in general suitable for nonlinear shallow- 
water waves; 

• We do not make any assumption on the velocity profile or on any related quan- 
tity. The only assumptions we make concern the value of the parameters e and 
fi (and a third parameter (3 linked to the amplitude of the bottom variations 
for uneven bottoms); the properties of the velocity profile (and in particular 
its vertical behavior) are then rigorously established; we thus believe that this 
article should clarify the discussions concerning the assumptions made on the 
velocity profile, and in particular the controversy raised in [9] about the Serre 
equations. 



2 General nondimensionalized water-wave equa- 
tions 

Parameterizing the free surface by z = ((t,X) (with X = (x,y) G R 2 ) and the 
bottom by z = —ho + b(X) (with ho > constant), one can use the incompressibility 
and irrotationality conditions to write the water-wave equations under Bernouilli's 
formulation, in terms of a potential velocity (i.e., the velocity field is given by 
v = V x , z # 

' d 2 + ^0 + (9 2 = O, -h + b<z<C, 

d n 4> = 0, z = -h + b, 

< <9 t C + VC-V0 = cU>, z = (, (5) 

k M + \(\ V0I 2 + (d z <P) 2 ) +g( = 0, z = (, 

where g is the gravitational acceleration, V = (d x , d y ) T and d n (p is the outward 
normal derivative at the boundary of the fluid domain. 

In deriving approximate equations by asymptotic methods it is necessary to intro- 
duce dimensionless quantities based on characteristic scales for the wave motion. 
Four main length scales are involved in this problem: ho the characteristic water 
depth, L the characteristic horizontal scale, a the order of free surface amplitude 
and B the order of bottom topography variation. Three independant dimensionless 



4 



parameters can be formed from these four scales. We choose: 



a 

h 



L 2 /X ' 



/l 



(6) 



where e is often called the nonlinearity parameter, while /i is the shallowness pa- 
rameter. 

Commonly two distinct nondimensionalizations are used in oceanography (e.g. Dinge- 
mans (1997)) depending on the value of \x: shallow water scaling [p, <C 1) and Stokes 
wave scaling for intermediate to deep water (ji 1). In this paper, we present a 
general nondimensionalization which applies to any wave regime. The order of mag- 
nitude of wave motion variables are given by the linear wave theory. In particular, 
(ghois) 1 / 2 and f^(^r) 1 ^ 2 are the characteristic scales of respectively the wave celerity 
and the potential velocity, with v = tant^/i 1 / 2 )//! 1 / 2 . 

Let us normalize all variables according to the scales anticipated on physical grounds: 



x = Lx', y = Ly' 

c = <, $ 



JL r / gho (jy 

ho 



z = h uz' , t 
b = Bb'. 



A 



Vghov 



t' 



From this general nondimensionalization we can recover the classical scalings for 
shallow and deep water: v ~ 1 when /i< 1 and v ~ /i -1 / 2 when [i ^> 1. 
The equations of motion ([5]) then become (after dropping the primes for the sake of 
clarity) : 



^ 2 /i<9 2 $ + + <9 2 $ = 0, 



'0, 



0. 



-u 2 fiV(-b) ■ v$ + 

V 

z/ 2 /iV(^C) • V$ + <9 Z $) =0 



d 4 $ 



' £ |V$| 2 + e 



-(-i + pb)<z<-t, 

^ = -(-l + /36), 
v 

v 



(7) 



2 V 



C = o, 



2 



C- 



In order to reduce this set of equations into a system of two evolution equations, we 
introduce the trace of the velocity potential at the free surface, namely ip = 3>| a=e ,„ c 
and the Dirichlet-Neumann operator as 



r£o, 2 _,e 



with $ solving the boundary value problem 



z/ 2 /i<9 2 $ + z/ 2 /i<9 2 $ + <9 2 $ = 



-(-1 < ^ < -C, 

V V 



(one can check that £/^[f(jV = \/l + |V(^C)| 2 ^n$|,_e c , where <9 n $ stands for the 
upwards nondimensionalized normal derivative at the surface). As remarked in 
[301ET], the equations (J7J) are equivalent to a set of two equations on the free surface 
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parameterization £ and the trace of the velocity potential at the surface ip = $\ z=e , v( 
involving the Dirichlet-Neumann operator. Namely, 

d t ( - ^g;[-c^ = o, 

ri,+r, ^ ivy/I- ^ (^[^ + ^V«)-V^) 2 ^ 
W + C + ^IWI 2(i + e vivci 2 ) = °- 

In order to derive the system ([9]), let Q ^[eC, , (3b\- be the Dirichlet-Neumann operator 
£/u[fC]' corresponding to the case v — 1. One will easily check that 

w>o, ^[ £C)/ 3fo] = I<^c], 

so that plugging this relation into ([8]) yields 

d t (-—g,le(,f3b]^j = 0, 

W + C + ^WI 2(1 + g2/i|vc|2) = 0- 

The system (|9]) recasts the water-wave equations in terms of the free surface elevation 
C and the velocity potential ip at the surface; this is the formulation which will serve 
as a basis for all the computations in this article. 



3 Shallow water models 

This section is devoted to the study of shallow water waves: fi < 1. The most 
general situation (when no assumption is made on e and (5) is addressed in §3.11 
where the the Green-Naghdi equations (also called Serre [32J or fully nonlinear 
Boussinesq equations [23]) are derived. The key point is an asymptotical expansion 
of the Dirichlet-Neumann operator {? M [eC] with respect to /i. 

We then briefly show how to recover simpler models under additional assumptions 
on £] the moderately nonlinear case e = O(^JJL) is addressed in §3.21 and the weakly 
nonlinear case e = O(fi) in §3.31 Some new models are then derived in the next 
sections. 

3.1 The fully nonlinear case: e ~ 1, j3 ~ 1 and 

From the assumption /i < 1, one gets v ~ 1 and we take i/ = 1 in (jHJ) to simplify 
(this corresponds to the usual shallow- water nondimensionalization). 

3.1.1 Asymptotic expansion of the Dirichlet-Neumann operator 

Since fi <C 1, we look for an asymptotic expansion of $ under the form 

N 

*«w ^/'''I',- (10) 

j=0 
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Plugging this expression into the nondimensionalized water-wave equations one can 
cancel the residual up to the order 0(fi N+l ) provided that 

Vj = 0, . . . , N, d 2 z <$> 3 = -d 2 ^ - d 2 ^ (11) 

(with the convention that $_i = 0), together with the boundary conditions 

V? = N 1 $ i U=.c = 5 o,i^, / 12 \ 
3 ' \ d z <& 3 = (3Vb • V^_i \ Z= _ 1+Pb {U) 

(where 5qj — 1 if j — and otherwise). 

Solving the ODE ffTTj) with (fT^J) is completely straightforward, and this procedure 
can be implemented on any symbolic computation software to compute the $j at 
any order. For our purposes here, we must take N = 1 and thus need to compute 
$ and $i; one finds 

$o = ip, (13) 
$i = (z-et)(-^(z + eC)-l + Pb)A,l> + l3{z-eC)Vb-Vil>. (14) 



Remark 3.1 This shows that at order 0(fi 2 ), and without assumption on e, the 
potential depends quadratically on z, as is well-known. When \x is not small, the 
expansion ( TJ7| ) is not convergent and this is the reason why we use another technique 
to expand the Dirichlet-Neumann operator in 

An approximation of order 0(fi N+1 ) of Q^[eC,, f3b]i( is then given by 

GvfeC, m = -^V(eC) • V$ QPP u=eC + d z <5> app lz=e( + 0^ N+2 )- (15) 
for N = 1, we thus obtain from (fTOl). ffT3]) and 

= -A*V • (W^) +/i 2 V • T[/i,/36] • + 0(/i 3 ) (16) 
where h = 1 + e( — (3b and where the linear operator T[h, (3b] is defined as 

T[h, (3b] -W = T*hT x ■ W + T*hT 2 ■ W. (17) 
with T* (j = 1, 2) denoting the adjoint of the operators Tj given by 

Ti ■ W := V • W - ^-(3Vb ■ W, and T 2 ■ W := -(3Vb ■ W. (18) 
V 3 2 2 

3.1.2 Derivation of the Green-Naghdi equations 

We use here the asymptotic expansion ffl6|) to derive the Green-Naghdi equations 
(also called Serre equations) from the full water waves equations (Q. 
One can equivalently write (fT6"j) under the form 

G^eCPbty = -/iV • (hv) + 0(/i 3 ), (19) 

where v = - n\T\h, (3b]Vip. 
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Remark 3.2 According to formulae (T3\)-(Lfi), the horizontal component of the ve- 
locity in the fluid domain is given by V(z) = V$o + V3>i(z) + 0(fi 2 ), and one can 
check that v = ^ Jf ^ +/3b V(z)dz + 0(/i 2 ), so £/ia£ ; m accordance with the classical 
derivations of the Green-Naghdi equations, v is the vertically averaged horizontal 
component of the velocity (up to 0(fi 2 ) terms). 

Recalling that v = 1 here, it follows from ( fl9l) that the first equation of ([9]) can be 
written + V • (hv) = 0(/i 2 ). 

From the definition of v given in (fl9l) . one gets easily 

W = v + /3b\v + 0(/i 2 ); 

taking the gradient of the second equation of replacing in the resulting 
equation by the above formula, using (I19j) . and neglecting the 0(/x 2 ) quantities 
gives the Green-Naghdi equations (recall that h = 1 + e( — (3b): 



d t ( + V • (/iv) = 0, 

h 



1 + ±T[h, Pb])dtv + VC + e(v ■ V)v (20) 



l*{ - 4 V ^ 3 (( V • V )( V ■ v) - (V • v) 2 )] + Q[h,(3b](v)} = 0, 



where 



T[h, pb}-W = T*hTi ■ W + T*hT 2 ■ W (21) 

= --V(/i 3 V ■ W) + (3-[V{h 2 Vb ■ W) - h 2 VbV ■ W] + (3 2 hVbVb ■ W 
3 2 

and the purely topographical term Q[h, 0b](v) (which is quadratic in v) is defined 

as 

Q[h,0b](v) = |-[V(/i 2 (vV) 2 6)-/ i 2 ((vV)(V-v)-(V-v) 2 )V6] 

+/3 2 ((v • V) 2 b)Vb. (22) 



Remark 3.3 The formulation [2D) of the Green-Naghdi equation is not at first sight 
the same as usual (see \2~$ and JEB$). A closer look shows however that they are 
exactly the same, as expected. The interest of the present form is that the second 
equation gives a straightforward control of the quantities (hv,~v), (/i7iv,7iv) and 
{h%y-, ^v) (with Tj given by ). This control yields regularizing effects of the 
same kind as those of the BBM equation compared to KdV, and is thus expected to 
ease numerical computations (work in progress). 

In [23], Wei et al. derived some Green-Naghdi equations (fully nonlinear Boussi- 
nesq models in that reference) with improved frequency dispersion by replacing the 
vertically averaged horizontal velocity v by the velocity v a taken at some intermedi- 
ate depth z a (x,y) (thus following the approach of Nwogu [15] for weakly nonlinear 
Boussinesq systems). Such systems could of course be derived similarly from (1201) . 
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3.2 The moderately nonlinear case: and e = 0{. s JJi) 

The simplifications that can be made on the Green-Naghdi equations f[2T)l) under the 
assumption e = O(yfp) depend on the topography. As for the surface variations, we 
distinguish three different regimes: 

(a) Fully nonlinear topography: (3 = 0(1). In this case, no significative simplification 
can be made, and the full equations must be kept. 

(b) Moderately nonlinear topography: (3 = 0( V //I). In this regime, the last term of 
the second equation of (12"0|) can be written 

-yV[((vV)(V-v)-(V-v) 2 )]+0(/i 2 ), 

so that one can replace ([20]) by 
d t C + V • (hv) = 0, 

(1 + ^T[h, (3b])d t v + VC + e(v ■ V)v - ^V[((v ■ V)(V • v) - (V • v) 2 )] = 0, 

(note that some simplifications could also be made in the term %T[h, /3b]dtv but 
they are not interesting because they would partially destroy the regularizing effects 
evoked in Remark 13. 3p . 

(c) Weakly nonlinear topography: f3 = 0(e). This stronger assumption allows a 
simplification of the term ^T[h, (3b]d t v, and one gets 

d t c + V • (hv) = 0, 

(1 - ^(^V-Dftv + VC + e(v ■ V)v - ^V[((v • V)(V • v) - (V ■ v) 2 )] = 0; 

the main interest of this model is that in the case of flat bottoms (b = 0), its 
unidirectional limit gives the Cammassa-Holm equation ([221 EH ESI ES])- 

3.3 The weakly nonlinear case: and e = O(fi) 

The assumption e = O(fi) is the classical long waves assumption which yields the 
usual Boussinesq models. Here again, we briefly show how to recover these models 
in different topographic regimes: 

(a) Fully nonlinear topography: (3 = 0(1). Neglecting in (1201) the terms which are 
of order 0(fi 2 ) under the assumption e = O(fi) yields the equations 

d t ( + V • (hv) = 0, 

(1 + f^T[h, (3b})d t v + VC + e(v • V)v = 0, 

equivalent to equations ([1]) and ([2]) and which correspond to the Boussinesq model 
for nonflat bottoms derived by Peregrine |12j . 

For many coastal applications the weak dispersion of these equations is a critical 
limitation. Several alternative equations, with the same approximations /i< 1 and 
e = 0(h), have been developed to improve dispersion properties (e.g. Madsen et al. 
P3], Nwogu p2], Schaffer and Madsen [16] or Bona et al. [T71H8]). 
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(b) Weakly nonlinear topography: [3 = O(ji). In this case, the equations simplify 
further into 

d t c + V • (hv) = 0, 

(1 - ^A)d t v + VC + e(v • V)v = 0, (23) 

for which topographic effects play a role only through the presence of h = l + eC,—(3b 
in the first equation. 

Remark 3.4 This latter Boussinesq system corresponds to the following approxi- 
mation of the Dirichlet-Neumann operator: 

S M K, W = ■ (>*W) - y v • AW + 0(/i 3 ), 

which is deduced from / TTPj) fry neglecting the terms which are 0(fi 3 ) in the present 
scaling. 



4 Fully dispersive models: e^/JI <^ 1 and e C 1 

Small steepness asymptotics (e^/JI <C 1) have been introduced by Matsuno [25] 
and generalized in [261 E3 EE] . It is often claimed that the derived equations are 
valid under the condition s^fji <C 1 only, and that the other asymptotic models 
(which satsify this condition) can be deduced from them. We show here that this is 
not always the case. The unified framework used in this article allows us to check 
for instance that the Matsuno equations do not degenerate into the Green-Naghdi 
equations in shallow-water and that their precision is much smaller. 
After making a new asymptotic expansion of the Dirichlet-Neumann operator in this 
physical regime, we also derive a new generalization of the Matsuno equations for 
3D flows with uneven bottom. 



4.1 Asymptotic expansion of the Dirichlet-Neumann oper- 
ator 

When e and (3 are small, it is possible to make a Taylor expansion of Q^[e(, (3b]ip 
with respect to e and (3\ such an expansion has been derived for one dimensional 
surfaces in [37]; this method could be generalized to two dimensional surfaces, but 
we chose here to use another technique based on the following formulas 

Kmi^eC^^-^fO,^) = -^[0,/36](C^[0,/3fo]^-//V-(C(VV)) 

Um^(^[0,^-^[0,0]V) = ^sech( v ^| J D|)[V-(6(sech( v ^I| J D|)V^))], 

where we used the Fourier multiplier notation: given two functions / and u, and 
denoting by^the Fourier transform, f(D)u is defined as: 

We only prove the second of these formulas because the first one can be established 
with the same techniques and can also be found in the literature (e.g. [38], and also 
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Th. 3.20 of [39] and Th. 3.1 of [29] where a similar formula is established for non 
necessarily flat surfaces). 

First recall that by definition, G^O, flb]i/j = d z <&[(3b}\ z=0 , where solves 

//A$[/%] + d%$\j3b] =0, in - 1 + (3b < z < 0, 
$[/36] u=0 = - M/5V$[/36] • V6), = ^ 

(we denoted instead of <£> to enhance the dependance on the bottom topog- 

raphy). We thus have, for all smooth function <p compactly supported in Q (with 
U={(X,z),-l + (3b(X) < z < 0}), 

(jjiA$[pb] + d%$[J3b])<pdXdz + f {$\J3b] - il>)<pdX 
a, Jz=o 



+ [ (d z ®[pb] - pfiVb ■ $ty3b])<pdX = 0. 

Jz=-l+3b 



Denoting x — lim^o ^(^[Z^*] — ^[0])> one can differentiate the above identity to 
obtain 

/ {fiA x + d\x)pdXdz + / x <pdX 
Jn Jz=o 

+ / (d zX -^b- V$[0] +6^ 2 $[0])v?rfX = 0. 



J z = -l 

Since this identity holds for all test function (p, and since <9^$[0] = — /xA$[0], we 
deduce that x solves the boundary value problem 

liAx + d\x = 0, in - 1 < 2: < 
XL=„ = 0, ^XU=- 1 = ^V-(6V$[0]|^_ 1 ). 

This problem can be solved explicitly: 

, N smh( x /uz\D\) V , 

< ( '- 2 '^^M' ,lw[1|| - Jl 

since moreover one has lim^o ^(^[0, f3b]ip—G^[0, 0]ip) = d z x\ z=0 and that 3>[0](-, z) = 
C ° S ^M%D\) Dl) ^> we § et the formula. 

A first order Taylor expansion of Q ^[e^ , [3b]ip with respect to e, together with the 
first formula, shows therefore that 

GpHiPW = S^Pbty - e^[0,/36](C^[0,/36]^) - efiV ■ (CVV) + O^e^ji) 2 )- 

(25) 

a first order Taylor expansion of £/ M [0, /?6]^> with respect to b, together with the 
second formula, gives also 

^[0, /3bty = G,[0, Oty + /i/?sech( v ^| J D|)[V ■ (6(sech(V^PI)VV))] + 0(V£(/V# 2 )- 

(26) 
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Let us now define the operator T M and as 

T M = - tanh y J|) V and B M = sech(^|D|)[&(sech(^|D|)-)] (27) 
We thus have £ M [0, O]-0 = • V^, and ([25]) and (J25J show that 

-e/iV ■ (CV^) + 0(^(e^Jl)\ v^(/3v^) 2 ). (28) 



Remark 4.1 In t/ie shallow water, weakly nonlinear regime, and with a weakly non- 
linear topography (that is, \i <C 1, £ = 0(h), [3 = 0(h)), one can deduce from K2B) 
that 

GM,PW = ■ W) - y V • AW + OGu 3 ), 
and we £/ras recover the approximation derived with another technique in § 3.3\ (see 



Remark\3.4\) 



4.2 Derivation of a fully dispersive model for 3D flows over 
uneven bottoms 

We derive here a new system of fully dispersive equations which generalizes to the 
case of 2D surfaces and nonflat bottoms the systems derived by Matsuno (ID sur- 
faces, flat [25] and uneven [26] bottoms), Choi [28] and Smith [10] (2D surfaces, flat 
bottoms). 

We first define the horizontal velocity at the surface as = (V$)|. =eC , where $ is 
the velocity potential given by fl24l) . By definition of ip and Q ^[e( , (3b]ip , we get 

W = V 5 + 5VC(<9 2 $)|_ C 

1 + e/j,|VC| 2 
= v 5 + £ y7i(T At -V5)VC + 0((£ v ^) 2 ), 

where we used ( )28|) and = vs + 0(e^fjl) for the last relation. Plugging this 
relation into ( 1281) . one gets similarly 



VP 



Q^eC, Pb\tl> = T M • v 5 + ■ (B^ s ) - e^JIT^ ■ (CVT M ■ v s ) 

-e^JlV • (Cv s ) + O^e^f, ^h(P^h) 2 ). (29) 

Taking the gradient of the second equation of and using the above two identities 
gives therefore the following set of deep-water equations: 



-TV • v 5 



d t v s + VC + ey/fl( 



yf\LV 



T M -(CVT M -v s ) + V-(Cv 5 )) 



(30) 



2JJlv 



v|v 5 r-vcT M -vc) = o, 
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where we recall that T M and are defined in (127)) . and that v = tanh( v //I)/ v //I (so 
that one can replace ^fjiv by 1 in f[3"Ul) in deep water). 

Remark 4.2 //we remove the topography term from these equations, we recover 
the two-dimensional equations (3.25)-(3.26) derived by Choi [28]. If we take the one 
dimensional version of (30}) we recover the equations derived by Matsuno (Eqs (19)- 
(20) of JMj for flat bottoms and (4.28)-(4.29) of (27] for nonflat bottoms). 

Remark 4.3 The equations [3D) are precise up to order 0(e+/Jl, E\fpi) in deep water 
(since v ~ n~ l l 2 ); one could also use them in shallow water, but they are then precise 
up to order 0(e,/3) only (since v ~ 1 in shallow water). This is the same accuracy 
as the one provided by the weakly nonlinear shallow water models (with a weakly 
nonlinear topography) of Section \3.3\ it is therefore not surprising to check that one 
recovers the Boussinesq system [23) from [3D) by a simple Taylor expansion of the 
operators and B^ and by observing that v = (1 + / u|A)v s + 0(fi 2 ). 
In the fully nonlinear regime (with fully nonlinear topography) studied in Section 
\3.1l the equations (3d}) are precise up to order 0(e,f3) = 0(1) while we saw that 
the Green-Naghdi equations $2~D) are precise up to order O(fi). It is therefore not 
surprising to check that the shallow water limit of (30}) does NOT give the Green- 
Naghdi equations [20}) . 

Remark 4.4 The system [3d}) is "fully dispersive" in the sense that its disper- 
sion relation is the same as for the full water-wave equations. This is the case 
because the expansion [29]) keeps the nonlocal effects of the Dirichlet-N eumann op- 
erator G^e^, /3b]ip . In fjlj \^2\ \43], the authors make a differential approximation of 
shallow water type based on Pade approximants; they show that the dispersive prop- 
erties remain good far beyond the shallow water regime when bathymetric changes 
are not too strong (in this latter case, the model is more complex and its range of 
validity much narrower [4$). 



13 



5 Conclusion 



In this paper we have presented a systematic derivation of the main 2D asymptotic 
models for shallow and deep water waves, which allows to clarify their validity 
domain. We have also derived a new 2D fully dispersive model, system (130]) . for 
small wave steepness which extends to uneven bottom the approach developed by 
Matsuno [27] and Choi [28]. We have shown that even though these models remain 
valid in shallow water, their precision is then far below what can be achieved with 
the Green-Naghdi or Serre model when fully nonlinear waves are considered (e ~ 1). 
Hence, contrary to what it is generally thought [251 E6], [27J [28] this approach cannot 
in practical be considered as a unified theory of nonlinear waves because it is not 
accurate enough for nonlinear shallow water waves (e~l, //<1). In particular, we 
have shown that system ( 130]) in the shallow water limit (/i< 1) does not correspond 
to the correct fully nonlinear equations, namely the Green-Naghdi or Serre equations. 
Another reason why these fully dispersive water models are not likely to furnish 
interesting models in shallow water is that there is no obvious shoreline boundary 
condition for them. 

The Green-Naghdi equations represent the appropriate model to describe nonlinear 
shallow water wave propagation and wave oscillations at the shoreline. For coastal 
applications, the Green-Naghdi equations can be easily extended to include accurate 
linear dispersive effects, which allow to describe shoaling processes in intermediate 
water depth (see Wei et al. [23] and Cienfuegos et al. [21], [5]). Another interest of 
these equations is that there is a natural shoreline boundary condition given by the 
flux conservation equation (the flux h\ vanishes at the shoreline in the first equation 
of [20]) . However, for deep water wave propagation the Green-Naghdi model is no 
more valid and the fully dispersive model (130]) must therefore be used. 
It follows from these considerations that the asymptotic description of coastal flows 
requires at least the use of 2 different models: one for shallow-water (e.g. the Green- 
Naghdi equations (T2TJT) ) and another one for deep water (e.g. (1311 ). The numerical 
coupling of these two models is therefore a natural perspective for further works. 
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